Calculations of the Local Density of States for some Simple Systems

A recently proposed convolution technique for the calculation of local density of states is described more thouroughly and new results of its application are presented. For separable systems the exposed method allows to construct the ldos for a higher dimensionality out of lower dimensional parts. Some practical and theoretical aspects of this approach are also discussed.Comment: 5 pages, 3 figure

Soliton effects in dangling-bond wires on Si(001)

Dangling bond wires on Si(001) are prototypical one dimensional wires, which are expected to show polaronic and solitonic effects. We present electronic structure calculations, using the tight binding model, of solitons in dangling-bond wires, and demonstrate that these defects are stable in even-length wires, although approximately 0.1 eV higher in energy than a perfect wire. We also note that in contrast to conjugated polymer systems, there are two types of soliton and that the type of soliton has strong effects on the energetics of the bandgap edges, with formation of intra-gap states between 0.1 eV and 0.2 eV from the band edges. These intra-gap states are localised on the atoms comprising the soliton.Comment: 6 pages, 3 figures, 3 tables, submitted to Phys. Rev.

Tight-binding modelling of the electronic band structure of layered superconducting perovskites

A detailed tight-binding analysis of the electron band structure of the CuO_2 plane of layered cuprates is performed within a sigma-band Hamiltonian including four orbitals - Cu3d_x^2-y^2, Cu4s, O2p_x, and O2p_y. Both the experimental and theoretical hints in favor of Fermi level located in a Cu or O band, respectively, are considered. For these two alternatives analytical expressions are obtained for the LCAO electron wave functions suitable for the treatment of electron superexchange. Simple formulae for the Fermi surface and electron dispersions are derived by applying the Loewdin down-fold procedure to set up the effective copper and oxygen Hamiltonians. They are used to fit the experimental ARUPS Fermi surface of Pb_0.42Bi_1.73Sr_1.94Ca_1.3Cu_1.92O_8+x and both the ARPES and LDA Fermi surface of Nd_2-xCe_xCuO_4-delta. The value of presenting the hopping amplitudes as surface integrals of ab initio atomic wave functions is demonstrated as well. The same approach is applied to the RuO_2 plane of the ruthenate Sr_2RuO_4. The LCAO Hamiltonians including the three in-plane pi-orbitals Ru4d_xy, O_a 2p_y, O_b 2p_x and the four transversal pi-orbitals Ru4d_zx, Ru4d_yz, O_a 2p_z, O_b 2p_z, are separately considered. It is shown that the equation for the constant energy curves and the Fermi contours has the same canonical form as the one for the layered cuprates.Comment: 21 pages, 10 figures, published in J. Phys.: Condens. Matter (complete and corrected References section

van der Waals interaction in nanotube bundles : consequences on vibrational modes

We have developed a pair-potential approach for the evaluation of van der Waals interaction between carbon nanotubes in bundles. Starting from a continuum model, we show that the intertube modes range from $5 cm^{-1}$ to $60 cm^{-1}$. Using a non-orthogonal tight-binding approximation for describing the covalent intra-tube bonding in addition, we confirme a slight chiral dependance of the breathing mode frequency and we found that this breathing mode frequency increase by $\sim$ 10 % if the nanotube lie inside a bundle as compared to the isolated tube.Comment: 5 pages, 2 figure

Dynamical properties of Au from tight-binding molecular-dynamics simulations

We studied the dynamical properties of Au using our previously developed tight-binding method. Phonon-dispersion and density-of-states curves at T=0 K were determined by computing the dynamical-matrix using a supercell approach. In addition, we performed molecular-dynamics simulations at various temperatures to obtain the temperature dependence of the lattice constant and of the atomic mean-square-displacement, as well as the phonon density-of-states and phonon-dispersion curves at finite temperature. We further tested the transferability of the model to different atomic environments by simulating liquid gold. Whenever possible we compared these results to experimental values.Comment: 7 pages, 9 encapsulated Postscript figures, submitted to Physical Review

Efficient index handling of multidimensional periodic boundary conditions

An efficient method is described to handle mesh indexes in multidimensional problems like numerical integration of partial differential equations, lattice model simulations, and determination of atomic neighbor lists. By creating an extended mesh, beyond the periodic unit cell, the stride in memory between equivalent pairs of mesh points is independent of their position within the cell. This allows to contract the mesh indexes of all dimensions into a single index, avoiding modulo and other implicit index operations.Comment: 2 pages, 0 figure

Core reconstruction in pseudopotential calculations

A new method is presented for obtaining all-electron results from a pseudopotential calculation. This is achieved by carrying out a localised calculation in the region of an atomic nucleus using the embedding potential method of Inglesfield [J.Phys. C {\bf 14}, 3795 (1981)]. In this method the core region is \emph{reconstructed}, and none of the simplifying approximations (such as spherical symmetry of the charge density/potential or frozen core electrons) that previous solutions to this problem have required are made. The embedding method requires an accurate real space Green function, and an analysis of the errors introduced in constructing this from a set of numerical eigenstates is given. Results are presented for an all-electron reconstruction of bulk aluminium, for both the charge density and the density of states.Comment: 14 pages, 5 figure

O(N) methods in electronic structure calculations

Linear scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high performance computers. The linear scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas is then discussed. The applications of linear scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear scaling methods are discussed.Comment: 85 pages, 15 figures, 488 references. Resubmitted to Rep. Prog. Phys (small changes